First-principles scheme for spectral adjustment in nanoscale transport 
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We implement a general method for correcting the low-bias transport properties of nanoscale systems within 
an ab-initio methodology based on linear combinations of atomic orbitals. We show how the typical problem of 
an underestimated HOMO-LUMO gap can be corrected, leading to quantitative and qualitative agreement with 
experiments. We show that an alternative method based on calculating the position of the relevant transport res- 
onances and fitting them to Lorentzians can significantly underestimate the conductance and does not accurately 
reproduce the electron transmission coefficient between resonances. We compare this simple method in an ideal 
system of a benzene molecule coupled to featureless leads to more sophisticated approaches such as GW and 
find a rather good agreement between both. We also present results for a benzene-dithiolate molecule between 
gold leads, where we study different coupling configurations for straight and tilted molecules, and show that 
this method yields the observed evolution of two-dimensional conductance histograms. We also explain the 
presence of low conductance zones in such histograms by taking into account different coupling configurations. 

PACS numbers: 73.63.-b,85.65.+h,72.90.+y,71.15.Ap 



INTRODUCTION 



Theories of coherent electron transport through nanostruc- 
tures are commonly based on a mean-field Hamiltonian H, 
describing a scattering region connected to crystalline cuiTent- 
caiTying leads. Starting from H, a scattering matrix S (or 
equivalent Greens function) is calculated and transport proper- 
ties such as the electrical conductance G obtained by evaluat- 
ing Landauer-type formulae. For example at zero temperature 
and bias, G = {26'^ /h)T{EF), where T{E) is the transmis- 
sion coefficient for electrons of energy E passing from one 
side of the scattering region to the other and E-p is the Fermi 
energy. The mathematical machinery for computing S and 
T{E) from H is well established and is often referred to as 
"the scattering approach" or "non-equilibrium Green's func- 
tion approach" to transport. When implemented with the same 
level of rigor, these are mathematically equivalent. Such the- 
ories have been used to predict transport propertiesof^a large 
number of molecules between metallic leads llll- 



The problem of identifying the most accurate mean-field 
Hamiltonian for predicting nanoscale transport is a more open 
question. Density functional theory (DFT) 1 10, 11] is a ver- 
satile tool for generating self-consistent, mean-field Hamil- 
tonians at both zero and finite applied voltages. For metal- 
lic systems, such as atomic chains between electrodes of the 
same material, there is quantitative agreement between ex- 
periment and predictions derived from DFT-based mean-field 
Hamiltonians. In contrast, for low-conductance systems such 
as single molecules attached to metallic electrodes, predic- 
tions can differ from experiments by orders of magnitude |5], 
due to limitations inherent to DFT |12-16]. These discrep- 
ancies mainly arise from an underestimation of the HOMO- 
LUMO (HL) gap, with theoretical predictions of the order 
of 40% of the experimental values lll7i [Tsll . Approximated 
exchange-correlation functionals also contain self-interaction 
errors lll4lll9[l . which are particularly important for localized 



states. Such errors, which worsen the agreement between 
the removal energy and the last occupied Kohn-Sham eigen- 
value, produce an incorrect alignment between the HOMO 
and the Fermi energy 1I20I1 and can also affect the width 
of the resonances when the coupling is small i24ll . Even 



in strongly-coupled molecules, such as molecules with thiol 
anchor groups, where the HOMO has a large weight on 
the sulphur atoms, localized states can still have large self- 
interaction errors and produce peaks very close to the Fermi 
energy. When the self-interaction correction (SIC) is ap- 
plied such states move downwards and the zero-bias conduc- 
tance decreases lf2lll . Other approaches that avoid the self- 
interaction error, at least partially, are hybrid functionals, such 
as the B3LYP 11221 12311 . that include part of exact exchange 
obtained from Hartee-Fock. Again, the main effect of such 
functionals on the transport properties is the opening of the 
HOMO-LUMO gap, specially when the molecule is strongly 
coupled to the leads and charge transfer is small [!24;]. Ex- 
act exchange approaches (EXX) such as Hartree-Fock or DFT 
with exact exchange [25] also open the gap ll24ll . A comple- 
mentary method to SIC is LDA+U ]j26, 27], where a parame- 
ter U is included to account for the intratomic repulsion, espe- 
cially that produced on the d states. The main effect of the U is 
again the opening of the HOMO-LUMO gap, which reduces 
the zero-bias conductance I28i1 . An alternative approach ap- 
plies Hedin's GW approximation on the extended molecule to 
coiTect its electronic properties ] 29 , 3^ 32,33]. GW is known 
to reproduce energy gaps more accurately and also takes into 
account screening effects of the leads. However, this method 
has technical difficulties related to the connection between the 
GW corrected extended molecule and the leads and is also 
very expensive. Consequently, to date, this method has only 
been applied to small molecules. 

Errors in the predicted spectrum of the DFT Hamiltoni- 
ans lead to erroneous predictions for the positions of reso- 
nances in T{E). To overcome this deficiency, spectral ad- 
justment in nanoscale transport (SAINT) can be used to ad- 
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just the peak positions of T{E) ^4, 35]. For a molecule 
attached to a metallic surface, the correct positions of the 
HOMO and LUMO levels could be determined experimen- 
tally from photoemission measurements. Alternatively, they 
can be estimated theoretically by noting that the HOMO (eh) 
and the LUMO (el) energies should coincide with the nega- 
tive of the ionization potential (IP) and electron affinity (EA), 
defined as IP = i?7v-i — En and EA = En — En+i, 
where A'^ is the number of electrons in the neutral molecule 



1 37|,|38t]. These quantities can be approximated by taking into 
account the Coulomb energy required to charge the molecule, 
Ec, as IP = -£h + Ec and EA = -el - Ec Q, so 
that the charging energy is equal to IP — EA = Ahl + 2i?c 
140], which is larger than the HL gap. However, when the 
molecule is placed in a polarizable environment such as metal- 
lic surfaces, the charging energy is reduced from the gas phase 
value by the polarization energies P_|_ and P_, which reduce 
both the HL gap and the Coulomb energy. The change AS 
of each frontier orbital can be decomposed into Coulomb- 
hole ( AEch), screened exchange (ASsx) and bare exchange 
or Fock (AEx) contributions 01611 . The bare exchange is in 
general small and the screened exchange satisfies ASsx ^ 
-2AScH for the HOMO and AEsx - for the LUMO, 
so that typically P-|_ ^ — P_ lfl6ll . From a physical point 
of view the reduction of the IP can be understood by noting 
that the work required to take one molecular electron to infin- 
ity is decreased, since the attraction of the positively charged 
molecule that the electron leaves behind is screened. The EA 
is increased because the positive image charges help stabilize 
the extra electron placed on the molecule. Finally, to account 
for all effects, in the presence of coupling to the leads, a fur- 
ther shift occurs due the coupling self-energy. 



THEORETICAL APPROACH 

Once the positions of the HOMO and LUMO resonances 
are known, the question arises of how to obtain the most ac- 
curate T{Ep). In the literature there are two approaches to 
this question and the aim of this paper is to compare the two. 
The first method selects either the HOMO or LUMO reso- 
nance e closest to the Fermi level and fits this to a Lorentzian 
T'{E) = V^/iiE - e)2 + r2], whose width F is obtained 
from the original DFT-based H. The zero-bias conductance 
is then given by G = {2c'^ /h)T' {Ep) Ijl^Q]. This method 
has been shown to give results in close agreement with exper- 
iments, but it is only reliable when the relevant orbital is close 
to the Fermi level and the transmission in the gap is described 
by the tail of a single Lorentzian. 

In our paper we compare this approach with an alternative 
method for achieving SAINT, which adjusts the diagonal el- 
ements of the mean-field Hamiltonian 134 14411 to reproduce 
known values of the IP and EA. The starting point is a pro- 
jection of the self-consistent mean-field Hamiltonian onto 
the atomic orbitals associated with the molecule, whose eigen- 
values and eigenvectors are {'i'n}n=i,...,N and {eri}n=i,...,Af, 



respectively, where N is the number of atomic orbitals on the 
molecule. From this it is possible to build a new Hamiltonian, 

= ffO + Ao ^ (*„J + Au ^ (*„J (1) 

no n^^ 

whose eigenvalues corresponding to the occupied and unoc- 
cupied levels are shifted by Aq and Au, respectively. Tak- 
ing into account the definition of the density matrix, p — 
Y.n=i I*") I fi'^n - Mc), where /ic is the chemical poten- 
tial, and the completeness relation, 1 = J2n=i 1^") '^^^ 
notes that in the limit of zero temperature /(ii^—/Xe) = 8(/io— 
E), and p = l*no) (*«J, so that E„„ {^nj = 

1 — /3. Finally, 



i? = H" + (Ao- A„)p + Aui 



(2) 



In the general framework of non-orthogonal basis sets, the ma- 
trix elements i/^^ of this Hamiltonian are obtained by adding 
to the matrix elements of H'^, the molecular matrix elements 
of the density matrix operator and the identity operator mul- 
tiplied by the respective constants. In the case of the density 
matrix one has 



N N 
n=l a,P = l 

where S is the overlap matrix. In this expression, the matrix 
elements p^p are defined to be pafi = Y.n=i CnQC*^/(e„ - 
/ic) where {c„q} are coefficients in the expansion \'^n) = 

The final result is a phenomenological Hamiltonian with 
matrix elements 



H^. = + (Ao - Au) (mI P W) + A„5^, 



(4) 



From this Hamiltonian one obtains a T{E) which possesses 
accurately-positioned resonances. In what follows, we com- 
pare the zero-bias conductance T{E-p) obtained from this 
method, with the value obtained by making a Lorentzian fit 
to the resonance closest to E-p. We also compare this method 
to most sophisticated approaches such as the GW approxima- 
tion, for a benzene molecule between featureless leads, and 
study how the qualitative and quantitative trends in the con- 
ductance are modified when this method is used to correct the 
transmission of a benzene-dithiolate molecule between gold 
leads. 



BENZENE RING BETWEEN FEATURELESS LEADS 

To implement the above SAINT, we utilize the quantum 
transport code SMEAGOL [45J, which uses the Hamiltonian 



provided by the ab-initio code SIESTA to compute self- 
consistently the density matrix and the transmission coeffi- 
cients. To illustrate how it works and compare it to other 
methods such as the GW method we study first the simple 
case of a benzene molecule coupled to featureless leads. We 
compare our results for the bias-dependent conductance to the 
results obtained by Thygesen and Rubio l32ll . 

In case of featureless leads, which act only as electron reser- 
voirs, the only effect of the coupling to the electrodes is a 
broadening of the molecular levels by the imaginary part of 
the self-energies, T. The retarded Green's function can there- 
fore be written as G^iE) = [{E + 2iT)S - H]-^. From 
the retarded Green's function we compute the lesser Green's 
function, both in and out of equilibrium and obtain the den- 
sity matrix to make the process self-consistent |45]. At the 
end of the self-consistent cycle, for a given voltage, we ap- 
ply the SAINT correction and calculate the density of states 
(DOS), the transmission, which in this case can be simplified 
toT{E) = (ir/2)Tr{[G'^(£;)-G«t s-}^ and the current. 
We define the SAINT shifts via the pai-ameter A = (Aq, Au). 

We used a double-(^ polarized basis set (DZP) for both car 
bon and hydrogen atoms and a real space grid defined with ai 
energy cutoff of 200 Ry. The DFT exchange and correlatioi 
energy was evaluated with the generalized gradient approxi 
mation (GGA) as parametrized by Perdew, Burke and Ernz 
erhof Ii47l] . With these approximations we obtain a HL ga| 
of 5.99 eV and a quasiparticle gap, calculated by energy dif 
ferences, of 10.00 eV, which is in very good agreement witl 
previous results |32]. 

We first analyze how to implement the SAINT method. Ii 
the case of a completely isolated molecule, the method is un 
ambiguous. However, when a molecule is coupled to elec 
trodes one has to separate the Hamiltonian and density matri? 
of the molecule from those of the electrodes [48]. The Hamil 
tonian of the whole extended molecule cannot be shifted, be 
cause that would introduce a discontinuity at the interfacs 
between the surfaces of the extended molecule and the bull 
leads. Therefore we project onto orbitals belonging to atom; 
of the molecule and define Hq to be the sub-matrix of the 
Hamiltonian involving these orbitals only. The molecular den- 
sity matrix is obtained by diagonalizing Hq. 

In Fig. ([T]l we compare results obtained by shifting the 
molecular Hamiltonian levels with those obtained by fitting to 
Lorentzians. Results are shown for the case A = (—1, l)eV 
and two different couplings to the leads, to illustrate again 
the effect of the broadening of the molecular levels. We fit 
Lorentzians to the nearest resonances and compare these to 
the T{E) obtained from the spectrally-adjusted Hamiltonian. 
As can be seen, simply by fitting to the nearest resonances, (ei- 
ther the HOMO or the LUMO), the resulting zero-bias trans- 
mission is much smaller than the SAINT result. Furthermore, 
adding both Lorentzians does not significantly improve the 
results. This disagreement arises because for a system with 
many resonances, the Breit Wigner formula does not accu- 
rately describe the electron transmission coefficient between 
resonances, even though it may be accurate in the vicinity of 




FIG. 1: Comparison between the full transmission curve corrected 
with a SAINT A = ( — 1, 1) eV and the transmission obtained by 
fitting the HOMO or the LUMO to Lorentzians (dashed hnes) and 
summing both (dashed-dotted hne). F = 0.12 eV in (a) and 0.25 eV 
in (b). Notice the different vertical scales on each panel. 




FIG. 2: Conductance obtained by numerical differentiation of the 
current (a) and the corresponding density of states at equilibrium (b) 
for various values of the SAINT correction. F = 0.25 eV. 



a resonance Oil] . 

We plot in Fig. (|2]) (a) the bias-dependent conductance ob- 
tained by differentiating the current. The uncorrected DFT 
result is very similar to the result of Thygesen and Rubio ll32ll . 
It has however slightly larger value, due to our use of more ba- 
sis states that produce additional resonances above the Fermi 
level, in contrast with the minimal and truncated Wannier 
function basis used in |32], which gives rise to only one reso- 
nance above the LUMO. We apply the SAINT correction with 
A = (-2, 2) eV and A = (-3, 3) eV, which roughly coiTe- 
spond to the correction to the HOMO and LUMO with and 
without taking into account the correction due to the image 
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charges, respectively. The former would also approximately 
correspond to the GW case and the latter to the Hartree- 
Fock case (HF). As can be seen, the SAINT corrected curves 
are significantly lower and flatter than the DFT result, which 
agrees with the GW and HF results. This change in the con- 
ductance can be explained by the opening of the gap and the 
subsequent reduction of the DOS in the middle of the gap, as 
can be seen in Fig. (|2]l (b), which reduces the transmission. 
Compared with 03211 . the absolute values are again slightly 
higher, which may be due to the use of different basis sets 
114 911 . Since the main effect of the GW and HF methods is 
the opening of the gap, which is exactly the same effect pro- 
duced by SAINT, the origin of the good agreement between 
these calculations is clear The GW and HF give rise however 
to additional structure in the occupied and unoccupied levels 
due to the energy-dependent GW self-energy and the differ- 
ent electronic structure method, respectively. Such structure 
can slightly modify the value of the transmission in the gap 
due to modifications in the tails of the resonances, but does 
not produce differences as large as those found between the 
exact result and the Lorentzian approximation. 

Two notes of caution should be added when the SAINT 
method is applied to more realistic systems. In some spe- 
cial cases, when the Fermi level is pinned at the HOMO or 
the LUMO and the charge on the molecule is significantly al- 
tered from its neutral value, the self-consistent charge transfer, 
corresponding to that given by the uncorrected DFT result, 
would introduce additional shifts in the levels and possibly 
additional structure. Such charge transfer would have been 
different if the gap had been bigger. Since the SAINT is non- 
self-consistent and the gap is opened at the end, the amount of 
charge transfer in such cases would be wrong. However the 
charge transfer in most of the systems is small ifsoll and large 
distortions due to it are not expected to be common. Another 
aspect that should be considered is related to molecular states 
located on the coupling atoms and which penetrate into the 
leads due to the strong coupling, such as the typical thiol-gold 
connections 1.341 . In these systems such molecular states have 
some weight inside the leads and therefore when the SAINT 
is applied to only the molecular orbitals, these states can be 
distorted. The main effect in the transport properties would 
be a change in the width of the associated resonances. These 
changes are not expected however to have an effect on the 
qualitative trends and would affect only slightly the quanti- 
tative predictions 115 ill . A special case, which would produce 
states that would penetrate deep into the leads, can occur when 
the coupling between the molecule and the leads is very strong 
and the on-site energies of both parts are very similar, but this 
is not very common in molecular electronics systems unless 
the molecule is very small 1,36.1 . 



BENZENE-DITHOLATE BETWEEN AU LEADS 

An archetypic example that has been extensively studied, 
both experimentally 0521 - 15811 and theoretically is the 




FIG. 3: A benzene-dithiolate molecule between (111) Au leads, tilted 
70° from the normal and contacted on the bottom surface in a hollow 
configuration and in the top surface to an adatom. 



benzene-dithiolate molecule coupled to gold leads. Recent ex- 
periments analyzed with 2D conductance histograms . 74 , 
75l l77ll . on similar molecules (oligo(phenylene ethynylene) 
(OPE) molecules capped with thiols) contacted on one side to 
a surface and on the other side to a scanning tunneling mi- 
croscope (STM), show a decrease of the conductance as the 
separation between the tip and the surface increases and a 2D 
zone with an upper and a lower edge of high and low con- 
ductances, respectively. Some histograms also show at large 
tip-surface separations a circular zone of low conductances 
142, TtIi . Such histograms suggest that these molecules can 
be contacted with a tilted configuration for small tip-surface 
separations; as the tip retracts from the surface, the angle in- 
creases and the conductance is reduced. This can happen, as 
we shall see, if at least one of the sulphurs is contacted to 
an adatom, which corresponds to the most probable coupling 
configuration to the tip. 

Many theoretical calculations have been carried out with 
the benzene molecule oriented normal to the surfaces and cou- 
pled in a hollow configuration, which is predicted to be the 
most stable contact configuration and gives a strong coupling. 
When the molecule is contacted on top of a gold atom and nor- 
mal to the surface, however, the coupling between the HOMO 
orbital, which is mainly made of perpendicular x and y p or- 
bitals, to the s orbital of the gold adatom is rather small due 
to symmetry. When the molecule is tilted, the x and y orbitals 
of the sulphur increase the overlap with the adatom orbitals 
and the coupling increases. The same does not happen how- 
ever when the sulphur is contacted in a hollow configuration 
because the coupling of such orbitals to the gold surround- 
ing atoms is already large and it is not severely reduced by 
symmetry. One would therefore expect a rather strong angu- 
lar dependence when one of the sulphurs is contacted to an 




FIG. 4: Transmission of a benzene-ditliiolate molecule between FIG. 5: Transmission of a benzene-dithiolate molecule between 
(ill) gold leads with both sulphur atoms contacted in the hollow (1 1 1) Soid leads with one sulphur atom contacted in the hollow con- 
configuration and the molecule normally-aligned (continuous line) figuration and the other on top of an adatom and the molecule posi- 
or tilted (dashed line). The results were obtained without (a) and ^1°"'=^ normally-aligned (contmuous Ime) or tilted (dashed Ime). The 
With (b) the SAINT corrections given in table l|l). 



adatom and a weaker angular dependence when both sulphurs 
are contacted in a hollow configuration. The additional low 
conductance zone in the histogram could be produced by an- 
other gold adatom on the surface, which would increase the 
distance between the tip and the surface and further decrease 
the conductance due to the small coupling on both sides. 

In order to test these assumptions, we carried out ab-initio 
calculations using the SMEAGOL code. The parameters were 
the same as for the benzene molecule, but in this case we had 
to explicitly include the leads. We chose (111) gold leads 
with 9 atoms per slice, 2 slices on each side of the extended 
molecule, and three additional slices to include the bulk leads. 
We also included periodic boundary conditions along the per- 
pendicular directions to make sure the transmission coeffi- 
cients were smooth. We chose a SZ basis set for the gold 
atoms, which included the s and d orbitals ll76ll . An example 
of one of the configurations, where the molecule is tilted and 
contacted to an adatom in one of the surfaces, can be seen in 
Fig. (O. We studied three coupling configurations, hollow- 
hollow (HH), adatom-hollow (AH) and adatom-adatom (AA), 
and for each case we compared the normally-aligned molecule 
(9 = 0° from the normal to the surface) to a tilted molecue 
(6' = 70°). 

Before showing the results it is necessary to comment a 
technical detail. In case of BDT and other molecules with 
degenerate levels, one has to take additional care in order 
to obtain the coiTect shifts. The HOMO resonance in the 
transmission coefficients given by this molecule is composed 
of two resonances, which correspond to two nearly degener- 
ated molecular orbitals, one of which is empty in the isolated 
molecule without leads [36]. These states become occupied 
when the molecule is saturated with hydrogens, which are 
added to the sulphurs, or when the molecule is coupled to the 
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FIG. 6: Transmission of a benzene-dithiolate molecule between 
(111) gold leads with both sulphur atoms contacted on top of an 
adatom and the molecule positioned normally-aligned (continuous 
line) or tilted (dashed line). The results were obtained without (a) 
and with (b) the SAINT corrections given in table l|l). 



leads. Consequently, if levels are shifted using the occupation 
obtained from the Hamiltonian of the isolated molecule, the 
empty level associated to the HOMO moves upwards, which 
is not the desired outcome. In such cases, to avoid this prob- 
lem, the LUMO has to be shifted downwards by Aq and only 
the LUMOh-1 and higher-energy eigenstates shifted upwards. 

In order to define the precise SAINT shifts it is necessary 
to know the DFT HOMO and LUMO, and the IP and EA (ob- 
tained by total energy differences) of the molecule in the gas 
phase, and the image charge corrections. The DFT HOMO, 
DFT LUMO, IP and EA are, respectively, -4.68, -1.38, 7.19 
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TABLE I: Distance between the image plane, located 1 A above the 
surface, and the center of the molecule (a), image charge correction 
(W) and final corrections to the occupied (Ao) and unoccupied (Ao) 
levels. 



TABLE II: Uncorrected and corrected zero-bias conductances G. 





a (A) 


WieV) 


Ao (eV) 


Au (eV) 


HH 0° 


4.34 


1.15 


-1.36 


1.45 


HH 70° 


2.68 


1.86 


-0.65 


0.74 


AHO° 


4.54 


1.10 


-1.41 


1.50 


AH 70° 


2.88 


1.73 


-0.78 


0.87 


AAO° 


4.74 


1.05 


-1.46 


1.55 


AA70° 


3.07 


1.62 


-0.89 


0.98 



and -1.22 eV, from where we obtain gas phase corrections of 
-2.51 and 2.60 eV, for the occupied and unoccupied levels, 
respectively. The precise gap varies however from config- 
uration to configuration, because it depends on the distance 
between the molecule and the image charge planes, which 
changes with the coupling configuration and the tilting angle 
of the molecule. We assume the image charge plane is 1 A 
above the surface and the point charge is in the middle of the 
molecule, so that the charging contribution can be approxi- 
mated hy W ^ (e2ln2)/(87reoa) fil, where a = d/2 - 1, 
and d is the distance between surfaces. The values of a, W 
and the final corrections to the occupied (Aq) and unoccupied 
(Au) states are shown in table 

The results for the HH configuration are shown in Fig. (|4]l. 
As said before, the HOMO peak in this system is composed 
of two almost degenerated resonances coming from bonding 
and antibonding molecular orbitals generated by the states of 
the sulphur atoms near the Fermi level I36ll . This peak is very 
high in energy due to self-interaction errors |f2^ and produces 
very high conductances, compared to experiments. When the 
molecule is tilted, the peak widens and moves down in energy 
due to the increase of the coupling |78]. Therefore, as can be 
seen in Fig. (|4]i (a), for 6 = 70°, due to the movement of the 
HOMO to lower energies, the zero-bias conductance of the 
normally-aligned molecule turns out to be higher than that of 
the tilted molecule. However, when the SAINT correction is 
applied, due to the facts that the transmission of the normally- 
aligned molecule decays faster than that of the tilted molecule 
and the shift of the occupied levels in the first case is bigger 
than the shift in the second case due to a smaller image charge 
correction, when the SAINT corrections are applied the situa- 
tion reverses and the normally-aligned molecule gives a lower 
conductance. The values of the transmission at the Fermi level 
can be seen in table dlll l. 

In the normally-aligned AH configuration, shown in Fig. 
(|5]l both the width and the height of the HOMO are smaller 
than those in the HH configuration. The decrease of the width 
of the resonance is due to the reduction of the coupling on the 
side of the molecule coupled to the adatom and the decrease of 
the height is due to the asymmetry of the contacts. When the 
molecule is tilted, the width of the resonance increases dra- 
matically and again moves down in energy. The LUMO is also 
affected this time by the increase of the coupling, which al- 





^uncorrected (Go) 


^corrected (Go) 


HH 0° 


3.3-10-1 


2.2-10-^ 


HH 70° 


1.310-1 


5.3-10-2 


AH 0° 


5.0-10-1 


4.5-10-^ 


AH 70° 


9.8-10-2 


2.4-10-2 


AA 0° 


8.1-10-1 


1.9-10-^ 


AA 70° 


1.2-10-1 


1.2-10-2 



lows a better communication between the states in the carbon 
rings and the states in the leads than in the normally-aligned 
configuration. At zero bias and without any correction, the 
conductance of the normally-aligned molecule is again larger 
than that of the tilted molecule, which does not agree with re- 
cent theoretical and experimental predictions 17711 . However, 
when the SAINT corrections are applied, the faster decay of 
the transmission and the larger shift of the occupied levels in 
the normally-aligned configuration again reverse the trend and 
the tilted configuration gives the largest conductance, as can 
be seen in table dllb . 

Finally, in the AA configuration, shown in Fig. (|6]l, the 
HOMO peak splits into a number of sharp resonances, one 
of which is pinned at the Fermi level. This movement of the 
HOMO towards the Fermi level in the AA or top-top (on top 
of surface gold atoms) configurations has been attributed to a 
reduction of the charge transfer, which moves the states up- 
wards and pins the HOMO at the Fermi energy Ii65i1 . This 
movement increases the zero bias conductance and, at the 
same time, decreases the transmission in the gap. When the 
molecule is tilted, the HOMO resonances merge into a single 
peak, whose width is much larger than that of the resonances 
of the normally-aligned molecule. This peak moves down- 
wards again and gives a lower transmission at the Fermi level 
than that of the normally-aligned case. However, when the 
SAINT corrections are applied we see again that the situation 
reverses and now the tilted molecule gives a higher conduc- 
tance. The results can be seen in table dUl i. 

The conductances observed in previous experiments range 
from 10-* Go S to 0.1 Go [56]. The values obtained by 
Reed et al. 1521] are rather small and could be produced by spe- 
cially stretched configurations. The rest of values are roughly 
within the range of conductances we obtain after applying the 
SAINT correction. The differences between the high and low 
conductance zones in different molecules observed in recent 
experiments range between approximately 1 order of magni- 
tude in pyridines 1421 and 2 orders of magnitude in OPE's 
1 7711 . In our case, which is similar to that of the OPE's, we can 
see a difference of roughly 1 order of magnitude between the 
highest and lowest corrected values (excluding the HH config- 
uration, which we predict would not be very probable in STM 
experiments). Probably the low conductance zone could have 
lower conductances arising from stretched AA configurations. 



7 



SUMMARY 

In summary, we have shown how to implement a spec- 
tral adjustment to the bare DFT-based mean-field Hamilto- 
nian, which corrects the positions of transmission resonances 
in T{E). We found that this method can improve both the 
qualitative and quantitative agreements with experiments and 
is an improvement over an alternative approach based on 
Lorentzian fits. 

For the case of a benzene molecule coupled to featureless 
leads, we tested the accuracy of using Lorentzian fits to cor- 
rect the transport properties of molecular junctions and found 
that such an approximation can severely underestimate the 
transmission in some cases. We also compared our method 
to previous GW results. We found a very good qualitative 
agreement, although the quantitative agreement was worsened 
by the presence of Lorentzian tails coming from unoccupied 
states. The agreement improved however when smaller num- 
bers of basis functions were used. 

We also used the SAINT method in a realistic junction of a 
benzene-dithiolate molecule between gold leads, in an attempt 
to explain the evolution of the conductance as a function of 
length in 2D histograms. We found that it was necessary to 
shift the molecular levels in order to obtain results that could 
explain the observed evolution, i.e. results where the zero-bias 
conductance of a normally-aligned molecule coupled to one 
or two adatoms was smaller than the zero-bias conductance 
of a tilted molecule. We also predicted the existence of a low 
conductance zone in the 2D histograms of these molecules lo- 
cated more than one order of magnitude below the main con- 
ductance zone. 

Finally, we remark that the SAINT method can be used 
to correcting mean-field Hamiltonians even when the IP, the 
EA and other corrections are not known from first principles, 
as is the case for long molecules such as those reported in 
J^]. To obtain a corrected Hamiltonian in these cases, the 
SAINT shifts Aq and Au can be treated as free parameters 
to adjust the transmission coefficient to yield agreement with 
some chosen property, such as the current-voltage character- 
istic, the temperature dependence of the conductance or the 
dependence of the conductance on a gate voltage in, for ex- 
ample, an electrochemical environment. Once Aq and Au are 
chosen to fit one of these measurements, the SAINT-corrected 
Hamiltonian can then be used to predict the others. 

V.M.G.S. thanks the Spanish Ministerio de Ciencia e Inno- 
vacion, the EPSRC and the Marie Curie European ITNs FUN- 
MOLS and NANOCTM for funding. 
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